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Abstract 

We present results on the full one-loop electroweak radiative corrections to the 
squark decay partial widths into charginos and neutralinos. We show the renormal- 
ization framework, and present numerical results for the third squark family. The 
corrections can reach values of ~ 10%, which are comparable to the radiative cor- 
rections from the strong sector of the model. Therefore they should be taken into 
account for the precise extraction of the SUSY parameters at future colliders. 



The Standard Model (SM) of the strong and electroweak interactions is the present 
paradigm of particle physics. Its validity has been tested to a level better than one per mil 
in the particle accelerators . Nevertheless, there are arguments against the SM being the 
fundamental model of particle interactions 0, giving rise to the investigation of competing 
alternative or extended models, which can be tested at high-energy colliders, such as the 
Large Hadron Collider (LHC) ||, or a e + e" Linear Collider (LC) jp. One of the most 
promising possibilities for physics beyond the SM is the incorporation of Supersymmetry 
(SUSY), which leads to a renormalizable field theory with precisely calculable predictions 
to be tested in future experiments. The simplest supersymmetric extension of the SM is 
the Minimal Supersymmetric Standard Model (MSSM) [|[. If the masses of the extra non- 
standard particles are very large as compared to the SM electroweak scale, the effects of 
these particles decouple from the SM low-energy effective Lagrangian |J. This means that 
if the extra particles are too heavy we could not discern between the SM and the MSSM 
by just looking at the low-energy end of the spectrum, since the only trace of the MSSM 
would be a light Higgs boson (M h o < 135 GeV) whose properties would not differ from 
the SM one. Nevertheless, when some of the extra particle masses are of the order of the 
electroweak scale, the next generation of colliders will be able to produce such kind of 
particles and investigate their properties. In this case non-decoupling effects appear || [| . 
While the LHC will be able to produce new particles with masses up to 2.5 TeV, the LC 
will be able to make precision measurements of their properties, provided they are not too 
heavy. For an adequate analysis, precise theoretical predictions are required, going beyond 
the Born approximation for SUSY processes. 

Up to now the major effort on the computation of SUSY radiative corrections has 
been put into the computation of virtual SUSY effects in observables that involve only SM 
external particles, or into the calculation of loop effects in the extended Higgs sector of the 
MSSMQ. But for the case of direct production of SUSY particles, one also needs a detailed 
knowledge of the higher-order effects for the processes with these SUSY particles in the 
external states. A number of studies have already addressed this issue, for production as 
well as for decay processes. For squark and gluino production in hadron collisions, the 
NLO QCD corrections are available |TT|; for squark-pair production in e + e~ collisions, the 



NLO QCD are also known, together with the Yukawa corrections |L2|- Concerning the 
subsequent squark decays into charginos/neutralinos, the QCD corrections were presented 
in || |13[0) whereas the Yukawa corrections were given in [[OJ. In this last work large 



corrections were found. They were derived, however, in the higgsino approximation for the 
chargino; hence, a full computation is required to consolidate the significance of the loop 
effects. 

In this letter we present, for the first time, a complete one-loop computation of the 
electroweak radiative corrections to the partial decay widths of squarks into quarks and 



1 See e.g. jL0[ and references therein. 

2 The gluino decay channel, which can be overwhelming for rriq > m q + nig, was studied in fjjfl . Here 
we will assume nig > g i. 
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charginos/neutralinos, 



(1) 



with q = b, t, q = b, t, x = Xa> Xi ■ We present the basic structure of the corrections, and 
illustrate their main features and their significance in representative numerical examples. 
The fine details, as well as a comprehensive analysis, will be presented elsewhere |J. Al- 
though the explicit results, for definiteness, are displayed for the third squark generation, 
our analytic results are valid for all kinds of squarks. 

In processes with exclusively SM particles in external states, it is possible to divide 
the 1-loop contributions into SM-like and non-SM-like subclasses. This separate treatment 
is often used in the literature, and it is useful since it allows to make the computation 
in small steps, checking each sector individually. As a distinctive feature of the radiative 
corrections to processes with supersymmetric particles in the external legs, this separability 
is lost. In this kind of processes the ultraviolet (UV) divergences of diagrams with virtual 
SM particles cancel the UV divergences of diagrams with non-SM particles. Any partial 
computation would yield UV-divergent and thus meaningless results. For this reason we 
have to compute all the SM and non-SM loop contributions at once, with the proper 
counter terms involving the renormalization of almost the full MSSM. 

We have used an on-shell renormalization scheme, which can be obtained by extending 
the on-shell scheme of jL6| to the MSSM. For the SM sector, the gauge-boson and fermion 
masses are treated as input parameters; the electron charge is defined in the Thomson 
limit, and the weak mixing angle is given in terms of sin 9yy = = 1 — M^/Mf. The SM 
sector and its renormalization is described in [I6[] and will not be repeated here. On the 
other hand, different conventions for the SUSY sector exist; therefore, a few comments on 
our procedure are in order. Concerning the Higgs sector, one mass has to be specified as an 
input quantity (for which we take the mass of the H ), and a definition for tan (3, the ratio 
between the vacuum expectation value of the two Higgs-boson doublets, t^/fi, is needed. 
Following [17], an indirect definition of tan f3 is given by the requirement^ SV2/V2 = 5vi/v\. 
For the sfermion sector, we use the input and the renormalization conditions as described 
in |fR)| ) which fixes the squark masses and the squark- mixing angles 9q. 

Finally, as a new ingredient, we have to specify the chargino-neutralino sector. The 
tree-level mass matrices are well known, but we list them in order to settle the conventions: 



M 



M° 



M 

V2M W cos (3 

I \r 




V2M W sin p 


M 



M z cos j3s w 
—M z cos (3cw 
Mz cos (3sw —Mz cos f3cw 
V — M z sin (3s w M z sin (3c w —fi 



-M z sin (3s w \ 
M z sin (3cw 







(2) 



M and M' are the SU{2)l and U{\)y soft-SUSY-breaking gaugino masses. These mass 



J For possible other definitions see e.g. 
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matrices are diagonalized by unitary matrices U, V, N via 



U*MV* = M D = diag (M l5 M 2 ) (0 < Mi < M 2 ) , 

N*M°N^ = M° D = diag (M?, M 2 °, M 3 °, M£ ) (0 < M? < M 2 ° < M 3 ° < M 4 °) . l J 

This sector contains six particle masses, but only three free parameters are available for 
an independent renormalization. As a consequence, we are not allowed to impose on- 
shell conditions for all the particle masses. For the independent input parameters, we 
choose: the masses of the two charginos and the mass of the lightest neutralino. By 
introducing counterterms for all the independent parameters in (Q), we are able to relate 
the counterterms of the fundamental parameters to the mass-counterterms SMi of the 
charginos. Similarly to we find: 



Mi SMx + M 2 SM 2 = M5M + fi5fi + SM^, 
M 1 M 2 (M 1 5M 2 + M 2 5M 1 ) = (M/x - M^ sin (2(3) 



MSfi + fiSM 



-Mir S sin(2/?) - sin(2/3) SM^ 



(4) 



The mass counterterms 5 Mi are fixed using the on-shell scheme relation, in the convention 
of [[HJ (but with opposite sign for E), 

SM, = -\ (tf L (M?) + Y? R (Mf)) - E^(Mf ) , (5) 

where E| Li?S | denote the one- loop unrenormalized left-, right-handed and scalar compo- 
nents of the self-energy for the zth-chargino. 8M' is determined from the lightest neutralino 
mass, inverting the relation 

NlJM%N* ip = 5Mi° , (6) 

where the neutralino- mass counterterm 5M® is fixed by the on-shell condition for Xi, i n 
analogy to (|5]) . It is a non-trivial check that with the counterterms determined in eqs. (|j) 
and (H), the one- loop masses for the residual neutralinos, computed as the pole masses, 
are UV-finite. The one-loop on-shell neutralino masses read 

Ml ° s = Ml + N* p 5M°^N* ai + i (E2(M Q 02 ) + E2(M° 2 )) + Eg(M° 2 ) , (7) 

where now the parameters of eq. (|2|) and the masses and mixing matrices computed in (^) 
have to be regarded as renormalized quantities. 

The choice of the lightest neutralino to fix the counterterm 5M' in is only efficient if 
it has a substantial bino component. If M' 3> (|/x|, M) then |AT n | <^ 1, and the extraction 
of SM' from (|(J) would amplify the radiative corrections artificially. In this case it would be 
better to extract SM' from the ath neutralino, such that |iVi Q | is large. This is, however, 
not relevant for the scenarios which are discussed in this letter. Notice also that our 
renormalization procedure makes use of positive-definite mass eigenvalues for charginos and 
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neutralinos, which require the introduction of some purely-imaginary non-zero elements 
in the iV-matrix (|j). Had we chosen a real iV-matrix, with some negative eigenvalues, 
the various renormalization conditions would be plagued with the explicit sign of the 
corresponding eigenvalue (see e.g. |20| ). 

At one- loop, also mixing self-energies between the different neutralinos and charginos 
are generated, which we write as follows: 

- l ±^(k 2 ) = -i (±f(k 2 )]tP L + ±f(k^P R + tf L (k 2 )P L + tf R {k 2 )P R ) ,a*i3, (8) 



with E denoting the renormalized two-point functions, and with the chirality projectors 
P{l,r} — 1 T 75- F° r the neutralinos, the renormalized self-energies are related to the 
unrenormalized ones according to 



n {L,R} — {L,R} ' ^SL 



Ef L -N ai 5M°; x N, x , tf R 



Xf R -N* ai SM° lX N; x , (9) 



and analogous expressions hold for the charginos. 

As far as vertex renormalization is concerned, the vertex counterterms are already de- 
termined by the renormalization procedure described above. In addition to the parameter 
renormalization, we have introduced also a field-renormalization constant for each left- and 
right-handed chargino and neutralino field. As an explicit example, we list the renormalized 
bottom-sbottom-neutralino vertex. The tree-level interaction Lagrangian reads |0| 1 



4e e kK{A^p L +A® a p R )b+h. c . 

V * a =l,2 »=!,.. .,4 



A+ a a — R\a (N* 2 — |%^al) — Africa N* 3 , — — \/2X b Ri^N a3 — ^ty^P^Nai , 

(10) 

where the bottom-quark Yukawa coupling is X b = mj/ (\/2M w cos (3). Introducing the one- 
loop counterterms analogously to |l5j we obtain the following counterterm Lagrangian Jl5 



9 z*,f,o 



d^yO hi 



^=hy° a (mJLp l + 5A m aa p R ) b + h.c. 



6a c 2 w (SM&r 8Mf 



+ sz a \ c x o abla + c x o bic sz c 



+ lv / 2 a a 
9 1* vf, 



AfJ- (5Z° Ra + 5Z b L ) P L + A® J- (5Z% + 5Z b R ) P R 



+ ^byl(6A^ aa P L + 6A^ aa P R ) b+ h.c. | , {c± 



5A 



(b) 



-V2 (5R { £\ b N a3 + 5X b R®N a3 ) + 5R { £ (iV a2 - h w N al ^j - Ut w R {b) a N al 



5A 



(6) 



-V2(SR^X b N* a3 + SX b R^N* 3 



'St 



^twN* x I —^-R^a + SR 2a 

6 \ tw 



(b) 



4 Note that our convention for the neutralino mass-matrices (0) is different from that of [Hi 



11 ! ""-^ ^ = z w (11) 

where 5a, <5M^ Z are the charge and mass counterterms for the MSSM, as given in pl| . 



The renormalized one-loop part of the amplitude for the decay b a — > 6%° can then be 
written as 



- £Z£* = [<W* + C_ a/3 P L ] v x o ; C ±a/3 = + Ai b L s + aS q 1pi . (12) 



Besides the counterterms 5 A from eq. (|TTD, it contains the one- loop contribution A 1PI to 
the one-particle-irreducible three-point vertex function and the quantity A s representing 
the contribution of the neutralino-mixing self-energies (H), 

A ±aa ~ 2^ A ±a(S M 02 - M 02 ' ^ ' 



Due to the presence of photon loops, the amplitude fll2|) is infra-red divergent; hence, 
bremsstrahlung of real photons has to be added to cancel this divergence. We therefore 
include in our results the radiative partial decay width T(b a — ► bXal)i including both the 
soft and the hard photon part. So finally, the complete one- loop electroweak correction is 
given by 

"o T o(b a ^bx° a ) u ° virt ^ r°(b a ^b x ° a ) ' 

/f m 2 _M 02 -m 2 VA Ci,) +A (b) C (b) )~2M°m,( A {b) C {b) +A (b) C (b) )\ 

0virt I (m 2 -Af02_ m 2 )( , A (f>) p+i^W |2)_ 4M m Re(A W ^W*) I ' 

for the neutralino decay channels, and by corresponding expressions 5°? for the chargino 
channels, as well as for the top-squark decays. 

The loop computation itself is rather tedious, since there is a huge number of diagrams 
to compute. This is better done by means of automatized tools. The computation of 
the loop diagrams has been done by using the Computer Algebra Systems FeynArts and 
FormCalc p2|,p3||. We have produced a set of Computer Algebra programs that compute 
the one-loop diagrams (and the bremsstrahlung corrections), which are then plugged into 
a Fortran code for the numerical evaluation with the help of the one-loop routines Loop- 



Tools p3fl . A number of checks has been made on the results. The UV and infra-red 
finiteness of the result, relying on the relations between the different sectors of the model, 
is a non-trivial check. We also have recovered results already available in the literature; 
for instance, we used our set of programs to reproduce the strong corrections of M, and, 



using the higgsino approximation, we could also reproduce the results of fL5| . Moreover 



we also checked that, when using the MS-scheme, the one- loop corrections to neutralino 



and chargino masses reproduce those of [19 



Although we consider the chargino and neutralino masses as input parameters, in our 
numerical study we treat them in a slightly different way. We choose a set of renormal- 
ized input parameters (M, M',/i), and apply (□), (131) to obtain the one- loop renormalized 
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Figure 1: (a) Tree-level branching ratios for the lightest sbottom branching ratios into 
charginos and neutralinos as a function of the lightest sbottom mass; (b) EW radiative 
corrections, eq. flT^|), to the partial widths of lightest sbottom decaying into charginos; (c) 
as in (b), but for the neutralino channels. 



masses. Of course, if SUSY would be discovered the procedure will be the other way 
around, that is, the MSSM parameters will be computed from the various observables 
measured, for example, from the chargino production cross-section and asymmetries at the 



LC p4j. For a consistent treatment, the one-loop expressions for these observables will 
have to be used |^ . 

As for the numerical analysis, we use a set of parameters relevant for the next generation 
of colliders. The squark masses have been chosen in the range of ~ 300 GeV. If the squarks 
have a mass around this scale, they will be produced at significant rates not only at the 
LHC, but also at a LC with 800 GeV center-of-mass energy. As for the gaugino sector, we 
make use of the GUT relation between the gaugino soft-SUSY-breaking mass parameters, 
M' = 5/3 tan 2 9wM. The input parameters for the SUSY-electroweak sector are chosen to 
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be 

tan p = 4 , n= -100 GeV , M = 150 GeV , M H ± = 120 GeV . (15) 
For the quark-squark sector we take 

m b = 5 GeV , m bi = 300 GeV , m b2 = m~ h + 5 GeV , % = 0.3 , 

m t = 175 GeV , m k = 300 GeV , 6~ t = 0.6 , (16) 

the rest of squarks are given a mass of ~ 1 TeV with a mixing angle 9q = 7r/4. All 
over our numerical results we apply the (approximate) condition for the non-existence of 
colour-breaking minima, demanding that the (computed) values of the trilinear soft-SUSY- 
breaking couplings A q do not exceed 3mq p6| . 

In Fig. |l] we present the results on the lightest-sbottom decay. Fig. |l]a presents the 
tree-level prediction for the various branching ratios {BR) as a function of the lightest- 
sbottom mass. We see that, aside from the third neutralino channel, all the channels have 
an appreciable branching ratio whenever they are possible. The opening of the bosonic 
decay channel {pi — > t\W~) is clearly visible at m bi = + My/ = 380 GeV. The 
charged Higgs boson channel (6i — > t\H~) opens at m~ b = mi x + M H ± = 420 GeV, but 
its partial width is much smaller (for the chosen value of the parameters), and hence we 
can not visualize its effect in Fig. |l]a. Whenever the bosonic decay channels are open they 
amount to a large fraction of the branching ratio [^7| , and they suffer from large radiative 
corrections |28| . In Fig. |I]b and c we present the one- loop electroweak radiative corrections 



to the chargino and neutralino channels, respectively. The lightest-neutralino channel is 
specially important, in that it is always open when we require the lightest neutralino to 
be the Lightest Supersymmetric Particle. When making a numerical scan over the MSSM 
parameters the corrections show a rich structure, owing to the large number of thresholds, 
pseudo-thresholds, etc. For example, the variation of the higgsino-mass parameter \i does 
not only change the chargino-neutralino masses and mixing angles, but also the value of 
the soft-SUSY-breaking trilinear squark couplings, and the mass of the heaviest top-squark. 
We see in Fig. [l]a the opening of the b\ — > channel at m~ b = 365 GeV. This threshold 
is accompanied by a corresponding divergence of the one-loop corrections to the lightest- 
chargino channel {5^ in Fig. |l|b), and also in the various neutralino channels (Fig. [TJc) . 
Of course, the (divergent) corrections near the threshold do not have a physical meaning, 
since perturbation theory breaks down. In the light-sbottom region {m b < 350 GeV) the 
branching ratio is neutralino-dominated. In this region the corrections amount to a 5-10% 
positive corrections for all of the channels. After the opening of the \2 channel the picture 
changes a little. The third neutralino can get corrections up to 30%, its BR, however, is 
smaller that 1%; these large corrections are thus of little interest. The rest of the channels 
continue with moderate corrections of the order of 5-10%; a special address to the heaviest- 
neutralino channel {xD with a 15% correction and a non-negligible BR ~ 1 — 8% is in 
order. At this point, however, we do not know yet the net effect on the corrected branching 
ratios since the electroweak corrections are similar in all decay channels. For quantitative 
statements also the QCD corrections H,|TB[ have to be included. 
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Figure 2: As in Fig. p] but for the heaviest sbottom. 



The corresponding tree-level branching ratios and radiative corrections for the heaviest 
bottom-squark are displayed in Fig. fj. They show a similar pattern to that of the lightest 
sbottom, but with more relevance for the chargino channels. We note that in this case neg- 
ative radiative corrections are attained for the lightest neutralino (Fig. ^|c). The maximum 
radiative corrections to the neutralino channel are 15% (for the heaviest neutralino), which 
has a non-negligible branching ratio all over the allowed sbottom-mass range. Notice that 
the very different corrections to the chargino channels (~ 17% for \i an d ~ 4% for X2) a ^ 
large sbottom masses will translate into corresponding corrections to the branching ratios. 

We now turn our attention to the lightest-stop decay channels. In Fig. [5] we present the 
branching ratios and the corrections to all possible decays as a function of the lightest top- 
squark mass. The chargino decay channels are important whenever one of the charginos is 
lighter than the top-squark. The top-squark is expected to be one of the lightest squarks 
of the model (due to large squark mixing), it could be even lighter than the top quark. In 
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Figure 3: (a) Tree-level branching ratios for the lightest stop branching ratios into 
charginos and neutralinos as a function of the lightest stop mass; (b) EW radiative correc- 
tions, eq. to the partial widths of lightest stop decaying into charginos and neutralinos. 



this latter case, if the chargino channel is closed, it would decay through Flavour Changing 
Neutral Currents, t\ —>■ c\\ [29[|- These two decay channels have been used for a experi- 



mental simulation to extract with high precision the top-squark parameters at the LC |p0| . 
For the chosen set of parameters the neutralino channels have never a BR larger that 20%. 
When the lightest top-squark is lighter than ~ 200 GeV the only possible decay channel 
has corrections to the width up to 20% (Fig. |3|b). These are, however, of little practical 
interest, since the BR is 100% in any case. Above this mass, both chargino channels re- 
ceive corrections in the 5-10% range. The neutralino channels obtain corrections in the 
5-15% range, but, as in the case of the sbottom decays, the largest corrections occur in the 
channels with smallest branching ratios. Again, several threshold structures are visible in 
the plot, accompanied by the corresponding divergence of the corrections. 

In summary, we have presented the set of full electroweak corrections to squark decays 
into quarks and charginos/neutralinos. These corrections can be sizeable, and therefore 
they have to be taken into account for the extraction of the MSSM parameters from 
experiment. A sample of the numerical results has been presented. The final impact on 
the branching ratios will also depend on the strong corrections [§,[13j . A combination of the 
two sets of corrections, as well as a comprehensive analysis will be presented elsewhere ||. 
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